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ABSTRACT 

In this note, I suggest that the model used to fit the X-ray surface brightness profiles of 
extended sources, like groups and clusters of galaxies, has to be corrected when the counts 
are collected in a wide energy band comparable to the mean temperature of the source and 
a significant gradient in the gas temperature is observed. I present a revised version of the 
/3— model for the X-ray brightness that applies to a intracluster gas with temperature and 
density related by a polytropic equation and extends the standard version that is strictly valid 
for an isothermal gas. Given a temperature gradient observed through an energy window of 
1-10 keV typical for the new generation of X-ray observatories, the f3 parameter can change 
systematically up to 20 per cent from the value obtained under isothermal assumption, i.e. 
by an amount larger that any statistical uncertainty obtained from the present data. Within 
the virial regions of typical clusters of galaxies, these systematic corrections affect the total 
gravitating mass estimate by 5-10 per cent, the gas mass by 10-30 per cent and the gas fraction 
value up to 50 per cent, when compared to the measurements obtained under the isothermal 
assumption. 
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1 INTRODUCTION 

The /3— model is widely used in the X-ray astronomy to parametrise 
the gas density profile in groups and clusters of galaxies fitting 
their surface brightness profile. Cavaliere & Fusco-Femiano (1976, 
1978) note that " ... since both gas and galaxies distributions con- 
forms to the same gravitational potential, the former can be directly 
related to the latter; the galaxies may be considered as tracers of the 
total potential well 
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Using the King approximation (1962) to the inner portions of an 
isothermal sphere (Lane-Emden equation in Binney & Tremaine 
1987; note that the King approximation is proportional to r^'^ at 
the outer radii, whereas the isothermal sphere is proportional to 
r~^cf. Figure 10): 

po,gai(l + a;^)"^''^, X ■ 



is the galaxies velocity dispersion and Tgas is the temperature of the 
gas. 

The surface brightness profile observed at the projected radius 
b, S{b), is the projection on the sky of the plasma emissivity, e(r): 

f°° e dr^ 

(Hereafter I adopt r as symbol for the projected radius b). 
The emissivity is equal to 

e(r) = A(rgas) "p erg s"^ cm"^, (5) 

where rip — pgas/(2.21/imp) is the proton density and the cooling 
function, A (Tgas), depends upon the mechanism of the emission 
and can be represented as 
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and the perfect gas law, one obtains the formula; 

Pgas = PO,gas(l + X ) ^' , 



(3) 



where the different energy distribution of the gas and galaxies is 
parameterized by using the parameter P ~ (crlj^i/Tgas), where cTgai 

* Throughout this note, a Hubble constant of 50 /150 km s~^ Mpc~^ is 
considered. 



where g is the velocity averaged Gaunt factor that is equal to 
about 1.2 within an accuracy of 20% for a bolometric emissiv- 
ity (for example, at Tgas > 2.5 keV, the emission is mainly due 
to bremsstrahlung and the cooling function can be written with 
A ~ 10"^^ and a = 0.5; see, e.g., Sarazin 1988). 

Assuming isothermality and a /3-model for the gas density 
(eq. ^), the surface brightness profile has an analytic solution; 
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Figure 1. The Binney & Tremaine (BT. from equation 4-125 in Binney 
& Tremaine 1987) dark matter profile for the self-gravitating isothermal 
sphere is here compared for different input parameters, [cr (km s^^), 
(or rs, Mpc)], to the Navarro-Frenk- White (NFW, 1995) profile that comes 
from extended and highly resolved numerical simulations of clusters of 
galaxies. Both of these are also compared with the King approximation 
(King 1962) to the inner part of the self-gravitating isothermal sphere. All 
of them are normalized to the central value of the self-gravitating isother- 
mal sphere profile [po = 9o-^/(47rGr^) = 9.05X 10~^^ g cm Inside 
the core radius, the NFW profile does not flatten like the BT profile. In the 
outer part of the region of interest (above 1 Mpc), agreement between the 
two profiles is obtained by increasing the velocity dispersion and the core 
radius (or scale radius) in the NFW profile. Fitting a power law, it can be 
shown that around 2.5 X rs ~ 2 Mpc the NFW profile approaches a r~^ * 
form. 



that is strictly valid under the condition that 3/3 > 0.5 and the 
coohng function A(rgas) does not change radially. 

In this note I will focus on energies typical for clusters of 
galaxies, considering the fact that there is now evidence for a de- 
crease in the gas temperature in the outer parts of clusters (Marke- 
vitch et al. 1998). Even if these results conflict with studies from 
other groups which indicate that clusters are generally isothermal 
(e.g. Irwin et al. 1999, Kikuchi et al. 1999, White 1999), I high- 
light how a temperature gradient can affect the estimate of the (3 
parameter. 

In the following discussion, I will assume that the cluster gas 
density, ngas, is well described by a (3— model and the gas is in the 
polytropic state, so that 

Tgas / rigas \ T ^ 

where the polytropic index 7 ranges between 1 and 5/3, the limits 
corresponding to the gas being isothermal and adiabatic, respec- 
tively. 



Figure 2. This plot shows the total dependence of the emissivity upon the 
temperature in different energy window [holometric, 1-10 keV, 0.5-2 keV] 
and for a cluster whose temperature ranges between 2 and 15 keV (where 
the bremsstrahlung emission dominates). The curves are calculated using a 
MEKAL model (Kaastra 1992, Liedahl et al. 1995) in XSPEC (version 10, 
Arnaud 1996) for two different values of metallicity: 0.3 (thickest symbols) 
and 1 .0 times the solar abundance. They are normalized to the bolometric 
emissivity at Tgas = 15 keV with Z = IZq. When the free-free ra- 
diation dominates, these curves can be calculated analytically integrating 
over the window energy [El — E2] the Gaunt factor, assumed equal to 
0.9(E/kT)-°-''^, multiplied by (fcT)-''-5e-(^''*^). Note that the emis- 
sivity is almost independent from the temperature of the plasma in the en- 
ergy accessible to the SOSAr observatory (cf. Table 1). 



Table 1. For a given energy window AE, I calculate the slope, a, of the 
power law in eqn. ^ over a selected temperature range ATgas . A is fixed 
to the value of the emissivity con'esponding to the higher Tgaa, i.e. a = 
log(e/e--'')/log{Tgas/Tg-J-). 



AE (keV) 


ATgas (keV) 


a{Z = 0.3Zq) 


a{Z = IZq) 


bol 


5-10 


0.45 


0.39 




5-7 


0.44 


0.37 




3-6 


0.43 


0.36 




3-5 


0.42 


0.34 


0.5-2 


5-10 


-0.13 


-0.20 




5-7 


-0.10 


-0.20 




3-6 


-0.08 


-0.21 




3-5 


-0.06 


-0.22 


1-10 


5-10 


0.25 


0.16 




5-7 


0.37 


0.27 




3-6 


0.47 


0.36 




3-5 


0.54 


0.41 



2 A POLYTROPIC /3-MODEL 

When the cluster temperature is above the energy range of the de- 
tector with a narrow bandpass (e.g. ROSAT), the emission measure 
will have a negligible dependence on the temperature (Figure 2). 

However, this dependence becomes significant when the en- 
ergy band is wide and its mean energy range is comparable to the 
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mean temperature of the cluster. For example, the emissivity due to 
free-free radiation for a plasma temperature of 10 keV is 50 per cent 
larger than one at 4 keV in the energy range [1-10] keV, whereas 
it changes by 1 per cent in the ROSAT bandpass [0.5-2] keV. In 
Figure 2, I show how the total cluster emissivity convolved with 
a given energy bandpass of a X-ray telescope depends upon the 
plasma temperature. 

A proper deprojection analysis (cf. White et al. 1997, Ettori & 
Fabian 1999) will be necessary to evaluate the emissivity in differ- 
ential volume shells and recover both the gas density and temper- 
ature profile. On the other hand, this technique is computationally 
very expensive and a simple fitting procedure can be preferred in 
most cases where the gas and total mass distributions are investi- 
gated. 

As discussed above, Cavaliere & Fusco-Femiano have intro- 
duced the model for an isothermal gas distribution. But if, as 
usually done, the gas density distribution is obtained from the pa- 
rameters of the best fit of eqn. ^ for the surface brightness profile, 
any temperature gradient present in the plasma will not be taken 
properly into account. To consider this correction, I assume that 
the gas density profile is well described by a /3— model (eqn. ^ and 
use the temperature dependence of the emissivity (given in eqn. 
^ through the parameter a) and a polytropic relation between gas 
density and temperature (represented by the index 7 in eqn. H) in 
the definition of the surface brightness (eqn. Then, eqn. ^can 
be re- written as 



^ Tg^s ?igas dr''^/ \/r'^ -b^ 
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2 r(3/j' -0.5) 2^0.5-3/3' 



5-0 



- 0.5-3/3(l+tii^) 



(9) 



Here /3' = P (l + a-^^) represents the uncorrected measured 
value. 

In Table 1, I calculate the a values for a set of interesting 
cases. 

It is worth noting that, even if the observed gas temperature, 
Tgas, is the projection on the sky of the real temperature, Tg^j' (^, 
weighted by the cluster emission, Tgas does not depend upon the 
parameter a (see also Markevitch et al. 1999): 
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(1 + ^2-j0.5-1.5/3[2+(7-l)(a+l)] 
(1 + x2)0-5-1-5/3[2+q(7-1)1 
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Figure 3. This plot .shows the corrections of the parameters fS (solid line) 
and 7 (dashed line) for a given temperature dependence of the emissivity, 
a, and an uncorrected polytropic index 7'. The thickest lines correspond to 
the case a = 0.6; the thinest lines are for a = —0.1. The other lines spans 
a values between and 0.5. 



(12) 



(13) 



/3' (7' - 1) = /3 (7 - 1) 
This system has the following solution: 

/3 = /?' (1 + a/2 - Q772) 

Y+a/2 — ci^'/2 

'y ^ l + a,/2-a7'/2 

Figure 3 shows the relative systematic corrections that affect 
the values of (/3, 7) for a given set of (a, 7'). 

In practice, once a temperature profile is measured in a known 
energy bandpass, the a parameter can be defined. The conversion 
from the count rate to the flux can be done assuming the central, 
highest temperature (i.e. with the largest emissivity, cf. Table 1). 
The functional forms, Sb = So{l + x'^y andTgas = ^o(l + a;2)^ 
can be then fitted to the surface brightness and temperature profiles, 
respectively. If they still represent a good model of the data, the 
correct (/3, 7) values can be estimated: 



1^ 0.5 — g — a6 

J = 1 + 



(14) 



The consequences of these corrections are discussed below. 



(10) 



I can then estimate the corrections that a polytropic tempera- 
ture profile produces on the uncorrected values (/3', 7'). To do this, 
I solve the system of equations given by the surface brightness and 
temperature profiles modeled with a /3— model: 



Soil + s2)«-5-3'3' = Soil + ^2)0-5-3,3(1+.^) 
To(l + 2.2 1.5,3' {7' -1) = + ^2)-1.5;3(7-l) 

that can be simplified to: 



(11) 



3 DISCUSSION AND CONCLUSIONS 

I suggest that the /3— model used to fit the X-ray surface brightness 
profiles of clusters of galaxies has to be corrected when the data 
from the next generation of X-ray observatories will be available. In 
fact, to include any temperature gradient that will affect the cluster 
emissivity as observed in a large energy window, we have to extend 
the use of the /3— model to the polytropic case. The new form of this 
polytropic /J— model is given in eqn. H. 
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Table 2. Changes of the derived quantities after that the corrections of 
the best-fit parameters are considered: unc(%) = 100 X (Qcon- — 

Qno — corr ) / Quo — corr ■ 



T 


A A J 




^7 gas 


a = 


0.47, 13 = 


2/3, 7 = 


1.20 


1 Mpc 


+8% 


-4% 


+12% 


2Mpc 


+13% 


-4% 


+18% 


''200 


+11% 


-5% 


+17% 


a = 


0.50, /3 = 


2/3, 7 = 


5/3 


1 Mpc 


+31% 


-10% 


+46% 


2 Mpc 


+56% 


-10% 


+74% 


''200 


+34% 


-10% 


+48% 



The spatially resolved surface brightness profiles obtained 
with the ROSAT observatory allow estimates of f3 with an accuracy 
(Icr) of about 2 per cent (e.g., Mohr et al. 1999, Neumann & Ar- 
naud 1999). However, the limited energy band pass of ROSAT does 
not arise any problem on the application of the /3— model when a 
temperature gradient is observed in the intracluster medium. This is 
also true for any surface brightness profile that is obtained collect- 
ing photons within an instrumental band pass at energies lower than 
the mean plasma temperature. In particular, considering that the ef- 
fective area of the present X-ray detectors is larger at £ ~ 1 — 2 
keV, an energy window around these values can be a good choice 
for a temperature-independent emissivity for hot clusters of galax- 
ies. But an energy range around few keV is still problematic for 
cool clusters and groups of galaxies that present temperature gradi- 
ents. 

The new generation of X-ray observatories, e.g., Chandrsj^and 
will operate in a wider energy band (e.g. [1-10] keV) than 
ROSAT and will provide a more accurate estimates of the param- 
eters of the /3— model. Therefore, as shown above, the presence of 
a plasma temperature gradient will affect the use of the model 
with a systematic uncertainty comparable or larger than any statisti- 
cal error. In particular, the estimate of (i) the total gravitating mass, 
Mtot oc /3 7 (1 + j;2)-i-5/3(T-i) (cf. eqn. (ii) the gas mass. 
A/gas oc ^(1 + x^)"^'^'' x^dx, and (iii) the consequent gas frac- 
tion, /gas, can be quite significantly affected. For example, given 
a cluster with a typical core radius, Vc, of 0.3 Mpc and a radial 
decrease of the plasma temperature from 6 to 3 keV with uncor- 
rected parameters (3' — 2/3 and 7' = 1.20 (e.g. Markevitch et al. 
1999), corrections of +8, -4, +12 per cent on Mgas, Mtot, /gas, 
respectively, will be necessary at r = 1 Mpc. At the more phys- 
ically meaningful radius r2oo, where the mean cluster density is 
A — 200 times the critical value (cf. eqn. ^ in Appendix), the 
corrections are +11,— 5, +17 per cent, respectively. These correc- 
tions increase considerably up to +56, —10, +74 per cent at r = 2 
Mpc for Mgas, Aftot, /gas, respectively, when a = 0.5, /3' = 2/3, 
7' = 5/3 (cf. Table 2). 
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APPENDIX A: OTHER ANALYTIC FORMULAE 

I write here the analytic expressions for the derived quantities like 
the total gravitating mass. A/tot, and the radius at which the mean 
density within a cluster at redhsift z is A times the background, 
TA, when the gas density is assumed well described by a (3— model 
and presents a polytropic (l<7<5/3) dependence upon the gas 
temperature: 



Mtot(r) 



r(r) r 



-( 



Gfim 
3 /37 To r, 



dlnp dlnT 

H 

9 In r 9 In r 



1.060 X lO^"* 



-P-yTorc 



(1 + 2:2)^ 



Mq, 



(Al) 



[Henriksen & Mushotzky 1986; see also Cowie, Henriksen & 
Mushotzky 1987, Hughes et al. 1988 for a critical discussion on the 
presence of an artificial cutoff in the parameter space, i.e. 2B < 3 
or 7 < 1 + 1/(3/3), to avoid the related virial density falling to 
zero (or negative values) at small radii] and, given that Mtot (''a) = 



|7rpc(l + ^)ViA, 
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3 13-7 To 



229.5 /37 To 
^i(l + 2)3 r|A 
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where the exponent i? = 1.5/3(7 — 1) + 1, To is the central tempera- 
ture in keV, rc the core radius in h^Q Mpc, fi is the mean molecular 
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weight in a.m.u. and the numerical values include the gravitational 
constant G, the mass of the proton rrip and all the unit conversions. 
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